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Abstract 

In this work we provide analytic results of infinite one-dimensional cellular automaton(CA). 
By realizing symbolic products, we investigate a subclass of infinite CA and prove analyti- 
cally that within this subclass the only allowed attractors are homogenous, steady and periodic 
states for arbitrary initial configuration. Our method also provide exact enumeration of these 
attractors and it is shown explicitly in a particular model. 

1 Introduction 

Cellular automata(CA) are simple mathematical models which can generate complex dy- 
namical phenomena. In principle, CA are discrete dynamical systems defined on a discrete lattice. 
The state of each site at any time t is in one of the g states. The interaction between sites is 
according to given local rules and the system evolves synchronously in discrete time steps. Cellular 
automaton was introduced by von Neumann in order to address the problems of self-reproduction 
and evolution ^J. However CA has not attracted much attention until John Conway introduced the 
game of lifes [2j around 1970 by using CA. About twenty years ago, Stephen Wolfram introduced 
cellular automaton to the physics community as models of complex dynamical systems (2i and a new 
approach to the parallel computing scheme(For a review see 0]). The interest in CA's potential 
application continues growing. In fact cellular automata have been used to simulate, for example, 
solar flares jSj, fluid dynamics j6j, crystal growth |7|, traffic flow and galaxy formation |^. 

Based on large number of numerical studies. Wolfram has suggested that CA can be classified 
into four classes. Cellular automata within each class has the same qualitative behavior. Starting 
from almost all initial conditions, trajectories of CA become concentrated onto attractors, the four 
classes can then be characterized by their attractors. According to Wolfram's classification, the 
class 1, 2 and 3 are roughly corresponding to the limit points, limit cycles and chaotic attractors 
in continuous dynamical systems respectively. More precisely their respective long time limits are 
: (1) spatially homogeneous state, (2) fixed (steady) or periodic structure and (3) chaotic pattern 
throughout space. The fourth class of CA behaves in a much more complicated manner and was 
conjectured by Wolfram as capable for performing universal computation. 

Due to the complexity of cellular automata, numerical computation of time evolution becomes 
the main scheme in CA studies. In P] some statistical approaches were introduced to characterize 
quantitatively the patterns generated by CA evolution. Statistical quantities such as entropies and 
dimensions can only provide average properties of cellular automata. However, exact results are 
always needed in more elaborate discussions . To make this point more explicit let us consider the 
effects of noise on CA. Due to the fact that all physical systems are coupled to noisy environment, 
it is interesting to know if the above classification of the attractors of CA remains intact under the 
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influence of noise. But this question can only be answered by knowing the deterministic CA exactly, 
and as a result the exact evaluation of the attractor without noise is called for. Unfortunately such 
exact calculations only exists for a few cases (A good review can be found in j^). For example, the 
GKL automaton was proposed in the 80' JT] and surprisingly the detail analytic proof was completed 
about a decade later For finite CA with majority rules, exact results have been obtained on 
finite lattice^Sj. It is shown in ^21 that the attractors are either steady states or periodic states of 
period two. It is important to note that the result of ^21 contains the majority rules as a special 
case. However, for a system of finite lattice, the configuration space is finite and can only contain 
periodic final states. Moreover, for infinite lattice, it is also known that the attractors are steady 
states if the initial state is finite|14|fThe finite initial condition will be defined in the next section.). 
Since the above results are only true for special configurations, the analysis for the attractors of 
infinite CA with arbitrary initial configuration is still lacking. It is the purpose of this work to fill 
this gap. 

In this work we provide the analytic results of infinite one-dimensional cellular automaton. We 
investigate a subclass of infinite CA and prove analytically that within this subclass the chaotic and 
complex structures are excluded. As a result, the only allowed attractors are homogenous, steady 
and periodic states. The main tool for proving the existence of steady states is the construction of 
a symbolic product among the fundamental blocks which will be introduced later. In section 2, we 
briefiy introduce notations and formulation of CA. In section 3 the product rules are given and 
several lemmas are proved to set the stage for the proof of the main result. The allowed attractors 
are shown explicitly for the case of fc = 2 in section 4. Finally, a brief discussion is given in section 
5. 

2 CA review and notations 

A cellular automaton is a spatial lattice of N sites where TV can be finite or infinite. The 
state of each site can be in one of g states at time t. Each site follows the same prescribed rules 
for up dating. For the rest of this work we will only concentrate on one-dimensional CA. The 
number of neighborhood of each site is denoted by 2k (In what follows k is always reserved specially 
for this context). The CA starts out with arbitrary initial configuration which is represented by 
^(0)— {*i5 ^2, S3, • • • , Sjv} where Si can be in any one of the g states. The configuration of system 
at time t is denoted by S(t)={si(t), 52(^)5 Ss{t), • • • , Sjv(t)}- In this paper the state of each site 
is restricted to g = 2 such that Sj 6 {1,-1} and the updating rule is given by the majority 
rule defined as follows. The neighbors of Si is defined as {si^a\a — ±1, ±2, ±3, • • • , ±k}. By 
introducing 

k 

a = — k 

the updating rule can be expressed as: 

r 1 if Qi >o 

Si(t + l)^l -1 <0 (2) 

[ Si{t) =0 

It is noted that the k — 1 case is also known as nearest neighbor CA. In this work, k is arbitrary 
and hence our results go beyond nearest neighbor CA. 

For finite one-dimensional CA, the lattice is arranged on a circle with periodic boundary con- 
ditions. Such cellular automata have a finite number of states, and as a result, after sufficient long 
time evolution the system must enter the state which can either be homogenous, steady or periodic 
state. Therefore the class 3 and 4 attractors do not exist in finite CA. In this paper exact results 
of a class of CA are established for infinite extend. 
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3 Lemmas and product rules 



To represent the state of each site more conveniently, from now on, the binary state {1, —1} 
wiU be graphically presented as {+, — } respectively. The one-dimensional CA starts out with 
some initial configuration which is represented by S(0)={si, S2, S3, • • • , Sjv} where each Sj can 
take either + or — . For any block of size m (with m > fc + 1) which contains only +(— ) site 
states is denoted as Am{Bm) (In the following discussion, the subscript of a block represents the 
block size, m is reserved to denote block size with m > fc + 1.): 

= {••• + + + •••} 

m 

= { •••---•• • }. 

m 

Given fc, any block of I sites without containing more than fc consecutive + or — site states is 
denoted by Xi. By using Am, Bm and Xi, any state can be symbolically decomposed as a product 
of Am, Bm and Xi, such as {• • • AmiXi^AmsBm^ • • • Xi^Bme • • •}■ For instance, for k — 2, 
the following configuration can be expressed as: 

{ + - + -- + - + + + -+} 

= { }{+- + -- + -}{+ + +}{-+} 

Obviously, this decomposition procedure is not unique. However, it is important to note that this 
decomposition procedure is just a tool for analyzing the state configuration at each time step and 
has nothing to do with the dynamics of system. During evolution, the size of each block may 
change and new block configurations will be generate. Our main idea of proving the existence of 
attractors is that the evolution of CA can he reduced to the evolution of these block variables. For the 
subclass of cellular automata considered in this paper, a fundamental properties which constitutes 
the foundation of the final proof are given as lemmas. 

Lemma 1 For 2k neighbors, the blocks A^^i and B^^i do not shrink during time evolution. 

The proof of this lemma can easily be seen by noting that the qi of each site in (Bfc_|_i) is 

positive(negative) or equal to zero. Hence these blocks can not shrink as the system evolves. In 
fact the size of any Am and Bm may grow during time evolution. This is due to the fact that all 
Xi configurations do not contain more than fc consecutive + or — site states, therefore the state of 
boundary site of Xi that is next to Am{Bm) will switch sign as the system evolves. In either case, 
the size of Am{Bm) grows. By assigning a notion of parity to Am as positive and Bm as negative, 
this growing process expands in both directions until both boundary sites meet another different 
parity block, then the size of the block remains the same in subsequent evolution. As a consequence, 
if the initial configuration contains either Am or Bm, then the structure of the attractor of these 
CA depends on the evolution of Xi. From the following lemma, any finite Xi can only lead to 
homogeneous or steady states, hence the chaotic state can never arise: 

Lemma 2 Any finite Xi with either Am or Bm o,s boundaries on both sides with m > fe + 1 will 
be eliminated during evolution. 

Since the lemma is quite obvious we omit its proof. In fact this lemma has been proved for any 
finite initial condition|14j which means that both Am and Bm are of infinite size. From Lemma 2 
one can see that these requirements are not necessary. Therefore, for any initial state without any 
block of Xao, the attractors of CA consist of the following configurations: homogeneous and steady 
states. For example with any fc, the state 
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»^ — -^ma ^7?^4 } 5 

is a steady state which does not change during time evolution. 

From above discussions, it is necessary to see if the chaotic state can arise from any initial state 
containing Xoo- By introducing and analyzing the combination of products of the fundamental 
blocks which will be defined below, the problem of can be reduced to a problem of finite Xi . 

For any Xi with I > 3(fe + 1), one can always pick a smaller section Xs(k+i) (C Xi) and 
consider the evolution of this block. This smaller section can then be split into three subblocks each 
of which has a size of fc + 1: 

^3(fc+i) — Xfc_|_iyfc_|_iZfc_|_i, 

For example, when k — 2, a, configuration X3(^k+i) — Xg — {— 1 1 1 — | } can be 

decomposed as: 

= {-- + - + -+ + -} = {- - +}{- + -}{+ + -}. 

In this way, any XQ(^f,_^_i^ block can be built up by all the possible Xfc_|_i's which are defined as 
fundamental blocks and we also denote them as X without the subscript. Obviously, by definition, 
the fundamental blocks do not contain more than k consecutive + or — site states. The number 
of fundamental blocks is 2(2'° — 1) and the total number of product combinations of fundamental 
blocks for constructing Xs(k+i) is bounded by 8(2'° — 1)^ which can be classified into five types: 
{XXX}, {YXX}, {XXY}, {XYX} and {XYZ}. Here, X, Y, Z represent different funda- 
mental blocks respectively. In order to proof the existence of the proposed attractors, one has to 
address how these possible configurations evolve in time. It is obvious that the number of product 
combinations grow rapidly as k increases. However, by realizing some relations among state config- 
urations, the number of product combinations needed to be considered can be reduced effectively. 

It is easy to check that the product combinations of {XXY}, {YXX} and {XXX} will 
produce Am or Bm which will not shrink in subsequent evolution. For example, for k — 2, the 

product of two fundamental blocks { 1 }{ 1 } will evolve to a configuration containing a 

block of B4 in the next time step. From Lemma 1 and 2, one can conclude that { — | — }{ — | — } 
can only contribute to the transient states. This situation will be addressed in detail for the case 
of A; = 2 in the next section. Therefore the cases required further investigation are the product 
combinations of {XYX} and {XYZ} types. 

As mentioned before, there exist relations among state configuration in one-dimensional CA such 
that the evolution of the product combinations can be simplified. The simplification is achieved by us- 
ing two types of operations which were pointed out in fJ . For any configuration S(t)={- • • s„_i, s„, 
• • •}, the conjugate state of S(t) is denoted as S{t) defined by flipping all the state of Sj. 

For example, the state S(0) = { h H 1 — } has a corresponding conjugate state S(0) = 

{H 1 Moreover, the state S(t) also has an image state associate to it. The image state 

S(t) of S(t) is defined by reversing sequence of the states of sites of S(t). That is to say, for any 
S(t) with N sites, S{t) — {si, S2, • • • , Sjv— 2; Sn— i: *jv}, the associate image state is 

S{t) — {Si, S2, • • • , Sn-2, Sn} 
— {^Nt Sat — 1: ^n — 2? ' ' ' t ^l}- 

For the S(0) considered above, the corresponding S(0) is { — | \--\ — }. It is important to realize 

that these operations are preserved under time evolution. These results are stated as theorems: 

Theorem 1 IfS{t) is the conjugate state o/S(t), then S(t + 1) is the conjugate state o/S(t + l). 

Proof: Under time evolution 

fc 



Si{t + 1) = sgn{ ^ (si+a + Si_„)}, 
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where sgn is the conventional sign-function and if J^i^i+a + Si— a) = 0, then Si{t + 1) = Si{t). 
The evolution of S(t) is given by: 



Si(t + 1) = sgn{ [si+a{t) + Si-a{t)]} 

a=l 
k 



k 



= Sgn{{-1) ^ [Si+ocit) + Si-acit)]} 
a=l 

When J2{si+a + Si-a) = 0, Si{t + 1) = Si{t), then 

Si{t+1) = i-l)Si{t) = Si{t). 

That is to say, if Sj(t) remains unchanged, so docs Si{t). The proof is complete. 

Theorem 2 //S(t) is i/ie image state o/S(t), i/ien S(t + 1) is the image state o/S(t + 1). 

Proof: The proof is similar to the theorem 1. The image of S{t) is given by 

^{t) = {S1,S2, • • • , Sjv-2, Sjv-1, Sn} 

with Si = Sjv— i- At the next time step, one has 

k 

Si{t + 1) = sgn{^[ 

k 

= sgn{ ^ [Sjv-(i+c.)(*) + Sjv-(i-a)(t)]} 



fe 



(*) + S(iv-i)-c«(*)]} 

= Sjv-t(* + l)- 

It is also true that if Sjv-i(i + 1) = Sjv-i(^), then Si{t + 1) = Si{t). Thus the image operation 
is persevered under time evolution. 

The method of establishing the existence of steady states is to show that all Xs(k+i) can produce 
Am or Bm as the system evolves. Once Am or Bm occurred, by using lemma 1 and 2, the steady 
state is obtained. From now on Am and Bm will be referred as the uniform blocks. However, 
one of the important properties of the conjugation operation is that the blocks Ama. and Bmb with 
ma = rUb form a conjugate pair. This fact is very useful in simplifying the analysis of the attractors. 

It is important to realize that the conjugation procedure divides the set of all possible product 
combinations of X^^k+i) into two mutually conjugate classes. Any X^^k+i) will belong to one 
of these two classes unless it is self-conjugate. The mutual conjugation between Am^ and Bm„ 
reduces the analysis on just one class. This reduction can be seen from the fact that if Q3(k+i) 
which is in either of classes produces a block Ama during evolution, then by applying theorem 1, 
the conjugate block Q3(fc+i) which belongs to the other conjugate-class will produce Bma.- As a 
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result, only one of the classes denoted by T requires further discussion. (The other conjugate- 
class is denoted by F.) On T, further reduction can also be obtained by considering the imaging 
operation. By considering two product combinations in T which form the image pair, if one of 
the pair produces either of the uniform blocks during evolution, then the other one results in a 
corresponding block respectively. Hence, one of these two product combinations can be eliminated 
without losing generality. However, in F, there are some product combinations that the image of 
them appear in T. Then the conjugation of these image product combinations are also in T and can 
also be eliminated. This fact can easily be seen by using the combined results of theorems 1 and 
2. The subset of T obtained after the above reductions is the set of product combinations which 
require detail analysis of the product rules. In the next section a detail discussion on the product 
rules for fe = 2 will be presented. 



4 The k = 2 Model 

For k = 2, there are six fundamental blocks: 

E = {+-+} G = {+ + -} H ={- + +} 

F = {- + -} / = {--+} J = {+--}, ^'') 
Prom (3), it is obvious that 

{E,G,H} (4) 

are the conjugate-blocks of {F, I, J} respectively. Moreover, the image operation produces the 
image pairs, (G, H) and (J, J). All Xg can be built up by the product combinations of (3) and 
altogether one has at most 216 product combinations. In considering the evolution of any Xg, as 
mentioned before, it can be checked explicitly that the product combinations of {XXX}, {XXY} 
and {YXXy types can produce the uniform block in the next time step and lead to the steady 
state. Note that this fact is also true for arbitrary k. Therefore, the main goal of the following 
discussion is to prove the existence attractors generated from {XYX} and {XYZ}. 

Let us start from the reduction procedures. It is noted that there are some products where Am 
or Bjn arises automatically as two fundamental blocks are grouped together. These products are 
excluded from the initial state with configuration X^o- For example, the product of E and G is 

{EG} = {+- + + +-} = {X2A3X1}. 

The complete list of these particular results is as follows: 

{EG, FI, GI, IG, HG, JI, HJ, JH, HE, JF}. (5) 

Moreover, there are also four products which can evolve to uniform blocks within one step and 

without referring to their neighboring fundamental block. These products are {EH, GE, FJ, IF}, 
and will also be omitted in the following discussion. As a result, the remaining products of two 
fundamental blocks are: 

{EF, EI, EJ, GF, GH, GJ, HE, HI, EE, EG, EH, IE, IJ, IH, JE, JG}. (6) 

This set can be divided into two subsets by conjugation relation. They are symbolically presented 
as follows for expressing the product relation more effectively: 
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(7) 
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F ^ 1 
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(8) 



Therefore, if the product combinations generated from (7) are in the class F, then the ones generated 
from (8) must belong to the corresponding conjugate-class T. Keeping these products and the 
conjugate relation between (7) and (8) in mind, the structure of attractors can be shown explicitly. 

Starting from {XYX}^ because of the conjugate relation between (7) and (8), one only need to 
use the group (7) to construct the product combinations of such type. By further considering the 
imaging operations, the product combinations can be reduced to only 4 cases: 

{EFE},{EIE},{GFG},{GJG}. (9) 

These combinations by itself do not contain enough information to justify that the system will 
produce uniform blocks at the next time step. However, it will be shown latter that these product 
combinations can lead to either the steady or periodic states. 

For the {XYZ} type, the analysis can be proceeded by concentrating on the ones in T. There 
are 14 product combinations required further consideration. By using the imaging operation for 
further reduction, the number of remaining product combinations is 10. Moreover, by combining 
the procedures of conjugation and imaging together, only 8 combinations remain: 

{EFG}, {EFH}, {EIJ}, {EIH}, {EJG}, {GFH}, {GHI}, {HFG}. (10) 

It is easy to check that except for the case of {GHI}, the rest of (10) will generate uniform block 
within one or two steps. For {GHI}, in order to check whether it can generate uniform blocks, it 
is necessary to take into account the effects of its neighboring blocks as the system evolves. Since 
{EG, IG, HG} is a subset of (5) which already contain the uniform blocks, the neighboring blocks 
of the leftmost side of {GHI} should not include any of {E,I, H}. Similarly, for the rightmost 
side of {GHI}, one should not include any of {F, G, J}. Thus, the following combinations require 
separate investigation: 

{FGHIE}, {FGHIH}, {FGHIJ}, {JGHIE}, {JGHIH}, {JGHIJ}. (11) 

By using the fact that the conjugate of {FGH} is {EIJ} which is already in (10), that is to say any 
block in (11) which contains {FGH} will generate the uniform block. This discussion also apply 
to {JGHIE} and {JGHIH}. However, the remaining case {JGHIJ} does not produce the 
uniform block. This particular block suggests a periodic configuration. In fact the infinite sequence 
of {JGHI} is a periodic state: 

{• • • JGHI JGHI ■ ■ .} 

i 

{• • . HI JGHI JG • • •} 

i 

{• • . JGHI JGHI • • •} 

Furthermore, by checking explicitly, any other neighboring block attached to {JGHIJ} always 
create a uniform block within two steps. 

Similar analysis can now be applied for analyzing the product combinations in (9). This can 
be illustrated by the case of {EFE}. According to the rules given by (7) and (8), the only 
neighboring blocks of {EFE} which do not generate uniform blocks is the infinite sequence: 
{• • • FEFEFEFE • • •} . This configuration is a steady state. By taking this block configura- 
tion with infinite sequence, the remaining product combinations in (9) are periodic states. For 
instance, considering {EIE}, the infinite sequence is {• • • EIEIEI • • •}. Under time evolution, 
this state oscillates: 
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{• • • EIEIEI • ■ •} 

i 

{■GFGFGF-} 

i 

{■ • • EIEIEI • ■ •} 

The other periodic state is: 

{• • • GJGJGJ ■ ■ ■} {■ ■ ■ IHIHIH ■ ■ ■} 

On the other hand, any other neighboring block attached to any blocks in (9) will produce uniform 
blocks by checking directly. In passing we note that the periodic attractors are period two as the 
finite CA. 

To summarize, for k — 2, from above analysis, we reduce the number of product combinations 
from 216 down to 12 ((9)+(10)) and it is shown that all attractors of the k — 2 model are consisted 
of three different types: 

1. The spatially homogeneous states such as A^o and Boo- 

2. The steady states of the form {...EFEF...} and {...AmiBmjAmj,Bmg...} with Vma > 
fe + 1. 

3. The periodic states which includes {...JGHIJGHI...}, {...EIEI...} and {...GJGJ...}. 

Therefore, the chaotic state and complex structure state do not exist in this cellular automaton. 

5 Discussions and conclusions 

By realizing fundamental blocks structure of the CA model with majority rule, it is shown in 
this paper that this model has certain nice properties which can be used to establish the structure 
of the attractors. It is shown explicitly that, with k — 2, for arbitrary initial configuration the 
only allowed attractors are the spatially homogeneous, steady and periodic states. Thus, the chaotic 
state and complex structure state referred in Wolfram's classification are excluded. Our method 
also provides explicit enumeration of all periodic states which are in fact period two as in the case 
of finite lattice ^21 ■ 

Although the proof was established for fc = 2, this approach is general enough for showing the 
result of arbitrary k. The approach is based on the fact that the uniform blocks Afc_(_i and Bk-\-i 
grow during evolution. Furthermore, by using the fact that any initial condition can be expressed as 
block products, therefore the problem of deriving the attractors turns into the problem of showing 
that any Xs(k+i) can generate the uniform blocks. Even though the number of product combina- 
tions of X3(^k+i) increases with fc, by using the conjugation and imaging procedures, the amount of 
explicit checking on the product analysis is greatly reduced. This approach has also been done for 
the case of fc = 3 but the analysis is too lengthy to be presented here and hence omitted. Thus, this 
work provides a complete analysis of the attractors of infinite CA with majority rules. Our results 
hold for any initial configuration in contrast to the results of previous works which assume either 
finite initial states or finite lattice. One can conclude that for this subclass of CA model, the only 
allowed attractors are homogenous, steady and periodic states. 
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